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Abstract 

In 1982 Hawking suggested that quantum fluctuations of the space-time metric will induce a non-unitary 
evolution from pure initial states pi to mixed final states pf in particle interactions at any energy. This 
hypothesis can be tested using existing CERN data on ir~p — > ir~ir + n on polarized target at 17.2 GeV/c. 
The purity of the final state pf is controlled by the purity of the recoil nucleon polarization. We develop 
spin formalism to calculate the expressions for recoil nucleon polarization for two specific measured initial 
pure states. Imposing the condition of purity on the recoil nucleon polarization we obtain conditions on the 
amplitudes which are violated by model independent amplitude analyses of the CERN data on polarized 
target at large momentum transfers. We conclude that pure states can evolve into mixed states in ir~p — 
Tr~Tr + n. In quantum theory such non-unitary evolution occurs in open quantum systems S interacting 
with a quantum environment E. The reduced density state is mixed and is given by Kraus representation 

Pf(S) = TrE(p/(S, E)) = ^2 SiPi(S)S^ . We show that Kraus representation leaves invariant the formalism 

1 

used in data analyses provided the co-evolution with the environment conserves P-parity and quantum 
numbers of the environment. The measured density matrix elements arc redefined to be environment- 
averaged density matrix elements. Measured elements that depend explicitely on the environment are 
predicted to violate certain phase relations. Measured elements that do not depend explicitely on the 
environment form a decoherence free subspace and satisfy the phase relations. This prediction is in excellent 
agreement with the CERN data and validates the view of pion creation processes as open quantum systems 
interacting with a quantum environment. The interaction can be thought of as a scattering of hadrons 
carrying energy-momentum and spin with particles of the environment carrying quantum entanglement. 
There is no exchange of energy-momentun between the hadrons and the environment in this non-dissipative, 
entanglement changing interaction. The observed process is time- irreversible and violates CPT symmetry. 
Following the proposal by Hawking, we suggest the origin of the environment is in quantum gravity rendered 
observable by the pion creation process. 
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I. INTRODUCTION: SPIN PHYSICS AND THE UNITARY S-MATRIX. 



The concept of scattering matrix was introduced by Wheeler in 1937 In 1943 Heisenberg in- 
troduced the concept of S-matrix to describe observables in interactions of elementary particles [3] . 
Heisenberg understood the S-matrix as a matrix of transition probability amplitudes between the 
initial and final states of the interactions. Just like the transitions of electrons between energy 
levels in atoms are not observable processes, the process of interactions of elementary paricles is 
not accessible to direct observation. Heisenberg believed that these classically unobservable pro- 
cesses are deeply related to our understanding of spacetime and will be eventually described by 
some sort of non-linear field theory that may encompass the irreversibility of quantum mechanics 
measurements 0]. In the absence of such a theory he proposed as a tentative solution the idea 
of asymptotic S-matrix defined in terms of in- and out- states in an analogy with potential scat- 
teringin quantum mechanics 0]. The unitarity of the S-matrix was demonstrated by Pauli in 
1946 |4| . Although Heisenberg believed that asymptotic S-matrix is a temporary tool, the success 
of Quantum Electrodynamics developed by Feynman, Schwinger and Tomonaga in 1940's provided 
evidence for its utility and meaningful physical interpretation. All Quantum Field Theories that 
followed were modeled on QED. To relate the theory to measurable S-matrix amplitudes, a re- 
duction formalism is first used to express the S-matrix elements in terms of product of vacuum 
expectation values of fully interactng fields 0, 0|. Next a unitary time evolution operator U(t) is 
found that takes the fully interacting fields to the free, asymptotic fields 

U(t) = U(t, -co) =Texp(- f dtHi(t)) =Texp(- f dt f d 3 xHi(x,t)) (1.1) 

J —oo J— oo J 

where T is time ordering operator and Hj(t) is the interaction Hamiltonian defined only with free, 
asymptotic fields. The S-matrix elements are then expressed in terms of vacuum expectation values 
of free fields which allows their calculation. The asymptotic S-matrix is defined as the limit 

S= lim U(t,-oo) = *7(+oo,-oo) (1.2) 

t^+oo 

The evolution operator U(t), and thus the unitary S-matrix, evolves pure isolated initial state 
into a pure isolated final state. Since the unitary evolution U(t) is time reversible, the interactions 
of particles must be time reversible isolated processes. 



In 1980 Wald showed rigorously that any scattering process of particles that evolves pure 
initial state into a mixed final state violates CPT symmetry and is time irreversible He 
suggested that such processes will occur in curved space-time, and that quantum gravity violates 
CPT symmetry and time-reversal invariance. In 1982 Hawking showed that in the presence of 
a black hole (macroscopic or microscopic) a pure initial state of interacting particles will evolve 
into a mixed final state as some of the quantum states produced in the particle interaction will 
fall behind the horizon and become inaccessible to quantum measurements by the observer jjj. 
He also suggested that quantum fluctuations of the space-time metric will have the same effect 
on interacting particles and induce their non-unitary evolution - at any energy 0, 0]. Hawking 
questioned a universal validity of the unitary time evolution (1.1) in the presence of metric 
fluctuations and suggested that initial and final state density matrices pi n and p ou t are connected 
by a linear but non-unitary evolution operator. Such mappings in fact exist and describe 
non-unitary evolution Pq,,±(S) = £{Pm(S)) of an open quantum system or process S interacting 



with an environment E 



11 



3 



Hawking's ideas inspired suggestions to test them experimentally. In 1984, Ellis, Hagelin, 
Nanopoulos and Srednicki proposed that quantum fluctuations of the metric form an environment 
with which interacting as well as free particles interact as open quantum systems 15J. Such 



interactions would lead to an observable loss of coherence and CPT violations in K°K° systems 
which maintain coherence over macroscopic distances. They supplemented Limblad time evolution 
equation for density matrix in dissipative open quantum systems with CPT violating terms 
to model interaction with quantum fluctuations and estimated the magnitude of their effects. 
Over the years and up to very recently other suggestions to test Hawking's ideas have been put 
forward 0, E3, E (3 M, EJ 0, . During the recent years experiments with neutral kaons have 
yielded remarkably sensitive results on violations of CPT symmetry and time reversal invariance, 
coherence of wave functions and entanglement of kaon pairs (24|. So far these experiments did 
not provide a conclusive confirmation of a non-unitary evolution of free kaon systems, and thus 
possible evidence for quantum gravity effects. 



In this work we return to the original idea of Hawking that pure initial state of interacting 
particles will evolve into a mixed final state as the result of interacctions with quantum fluctuaions 
of the metric during the process of particle interaction. To test unitarity in particle scattering 
requires to examine a variety of exclusive hadron processes with known initial and final state 
density matrices measured in spin physics experiments using polarized targets or polarized beams. 
Spin physics experiments make no use of the unitarity assuptions (1.1) and (1.2) and can thus 
inform us about their validity. 



Spin physics as a research field was initiated in 1946 and 1948 by Julian Schwinger. Schwinger 
wanted to know if nuclear forces are indeed spin independent as was generally assumed at that 
time. He designed experimental methods to create polarized neutron beams and to measure their 
polarization 25|, |26J]. In 1949 Wolfenstein introduced the concept of spin observables measurable 
in Schwinger 's experiments with nucleons 27[. In 1957 Chamberlain and collaborators at Berkeley 
reported the first evidence for spin dependence of proton-proton and proton-neutron scattering at 



320 MeV 



An important element of their experiments was the measurement of recoil nucleon 



polarization in a secondary rescattering, a method first proposed by Schwinger. 



In 1961 Bethe and Schumacher introduced the concept of amplitude analysis - the construction 
of scattering matrix from the measurements of a complete set of spin observables [29]. In 
1972 van Rossum, an earlier associate of Chamberlain at Berkeley, and his spin physics group 
at CEN-Saclay reported the first measurements of recoil nucleon polarization in ir + p elastic 
scattering at 6 GeV/c at CERN [30]. These measurements closed the set of spin observables 
and enabled them to perform the first amp litude analysis of a hadronic reaction. The measured 
pion-nucleon amplitudes at 6 GeV/c [3l( invalidated all Regge models, and emphasized the 
importance of experimental measurements of scattering amplitudes in experiments using polarized 
targets. In 1978 Lutz and Rybicki extended the concept of amplitude analysis to pion production 
processes ttN — > tt~it + N and showed that a nearly complete amplitude analysis is possible using 
only measurements on transversely polarized target 32J. Their results were later generalized to 

33L 34 ] . Amplitude analysis formalism was also extended 
36J]. A comprehensive introduction to spin in particle physics and 
modern spin physics technology is given in Ref. 37( and in the recent book by Leader [38j. 



7r p — > 7r°7r°n and t p-> irmi >i:orossos 
to inclusive measurements ]3 



The result of any possible experiment of a considered reaction \i >^ \f > involving measure- 
ments with spin of particles is the measurement of its final state density matrix pf. Given the 
initial state density matrix pi, the final state pf is given by a quantum evolution of the initial state 
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Pi determined by the 5-matrix 



PS = S Pi S + (1.3) 

The equation (1.3) allows us to express pf in terms of the full set of transition amplitudes 
describing the reaction. At first we are interested in two processes n~p — > ir~p and ir~p — ► 7r~7r + ra. 
As we shall see later, the purity of the initial state Pi{it~p) is controlled by target polarization 
vector P while the purity of the final states pf{ir~p) and pf(ir~TT + n) is controlled by recoil nucleon 
polarization vectors Q{P) and Q(Q, P), respectively. Here = (0,<p) specifies the direction of 7r _ 
in the center of mass systems of the two pions. The initial state is pure if and only if \P\ 2 = 1. 
Similarly, the final states are pure if and only if the recoil polarizations have a unit magnitude 
|Q(P)| 2 = 1 and \Q{fl,P)\ 2 = 1 0- 



Using known expressions for recoil nucleon polarization in n~p — * n~p 39|] we find that pure 
states evolve always into pure states. The same conclusion holds true for a class of meson baryon 
two-body reactions with the same spin structure. Such is not the case for ir~p — ► tt~7t + ti. 



Following the initial work of Lutz and Rybicki [32], we develop a spin formalism in 
ir~p — ► ir~ir + n to express the final state density matrix and recoil nucleon polarization in terms 
nucleon transversity amplitudes with a definite dipion spin, helicity and naturality. We then 
calculate exact expressions for |Q(S7,P)| 2 for two special pure initial states with transverse 
polarizations P y = ±1 and impose the condition that the corresponding final states is pure. For 
each P y = ±1 we obtain an elegant relationship between moduli and certain relative phases of 
all amplitudes that must be both satisfied by the unitary S'-matrix. These unitarity conditions 
are satisfied for dimeson masses where only S- and P-wave contribute. At higher masses the 
unitarity conditions require that either (A) relative phases between natural (unnatural) exchange 
amplitudes of the same nucleon transversity are all zero or (B) natural (unnatural) exchange 
amplitudes with dimeson spins larger than 1 all vanish. These conditions must hold true at all 
energies s, dimeson masses m and momentum transfers t. The unitarity conditions (A) and 
(B) also apply to a class of similar processes with dimeson final states that have the same spin 
structure, such as ttN — >■ 7r°7r°iV, KN — >■ KttN, ttN — >■ Kit A etc.. 

CERN measurements of tr~p — > 7r~7r + n on polarized target at 17.2 GeV/c find non-zero 
unnatural exchange amplitudes in S-, P- and -D-waves and non-zero natural exchange amplitudes 



in P- and D-waves at small 42], |43J as well as at large 44j] momentum transfers t. These 
measurements also find non-zero phases between all unnatural exchange amplitudes and between 
the P-wave and D-wave natural exchange amplitudes, at small and large t. Both unitarity 
conditions (A) and (B) are thus violated by the CERN measurements. We conclude that in 
7r~p — * n~ii + n and in similar pion creation processes pure states can evolve into mixed states in a 
non- unitary evolution of initial states pi into final states pf. 

In quantum theory non-unitary evolutions occur in open quantum systems interacting with an 
environment 



111 . 112 . Il3l . Il4| . The co-evolution of the system S with a quantum environment E is 
a quantum operation pj(S,E) = U Pi(S, E)U + with a unitary evolution operator U. The initial 
state of the combined system is assumed to be separable state pt(S,E) = pi(S) (g> Pi(E) but the 
the final state pf(S,E) is an entangled state of S and E. The measured state pj{S) is a reduced 
state given by Kraus represetation 

p f (S) = Tr E ( Pf (S, E) = J2 S iPi (S)S+ (1.4) 
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where Sg =< eg\U\eo > are operators acting on the state space of the system S, \e% > are 
quantum states of the environment, and where we assumed that the initial state is a pure state 
Pi(E) = | eo >< eo | - Kraus representation describes a non-unitary evolution from Pi(S) to p/(S), 
and the state pj(S) is a mixed state. 

At first sight there appears to be a contradiction. Starting with the unitary evolution (1.3), a 
spin formalism is developed to analyze the measured data. The data analysis leads to conclusion 
that the evolution is non-unitary and described by Kraus representation (1.4). We show that 
Kraus representation leaves the form of all equations invariant provided that the co-evolution 
U conserves P-parity and quantum numbers of the environment. However, the measured 
density matrix elements and the measured moduli and phases are environment-averaged values 
of co-evolution density matrix elements and co-evolution amplitudes that explicitely depend 
on interacting degrees of the environment. The conclusion, that the measured density matrix 
elements are averaged values predicts a violation of certain phase relations. We show that such 
phase relations are violated by the CERN data on ir~p — > n~ir + n. 

The CERN data validate the Kraus representation and the new view of pion creation processes 
as open quantum systems interacting with a quantum environment. We can think of the 
co-evolution process as a scattering of initial hadron states that carry energy-momentum and spin 
with particles of the environment that carry quantum entanglement. There is no exchange of 
energy-momentum between hadrons and the environment in this non-dissipative process. Since 
environment states do not have well defined antiparticle states and the interaction is non-local, 
the coevolution process violates CPT symmetry. Since the measured process is non- unitary, 
it is time- irreversible and violates CPT symmetry. Following the proposals by Wald [7j and 
Hawking [1, 0] , we suggest that the observed quantum environment originates in quantum gravity. 
In this view, pion creation processes act as non-classical instruments that make observations of 
quantum gravity possible. 

The paper is organized as follows. In Section II. we present the form of the final states and 
conditions for its purity in rr~p —* 7r~7r + n and develop the necessary spin formalism that will be 
used in this work and in sequel papers. In Section III. we show that pure states evolve into pure 
states in ttN — > ttN and in similar processes. In Section IV. we derive conditions that must be 
satisfied by the amplitudes in n~p — > TT~ir + n in order for two specific pure initial states to evolve 
into pure final states. In Section V. we show that these conditions are violated by CERN data on 
ir~p — > TT~TT + n on polarized target and conclude that in ir~p — > 7r~7r + n pure states can evolve into 
mixed states. In Section VI. we comment on a formal connection between a change from unitary to 
non-unitary evolution and a change in metric. In Section VII. we briefly review the basic concepts 
of evolution of open quantum systems interacting with environment and Kraus representation. In 
Section VIII. we show that Kraus representation leaves invariant the spin formalism developed 
in Section II. provided that the co-evolution with environment conserves P-parity and quantum 
numbers of the environment states. The prediction that the measured density matrix elements 
and bilinear terms of amplitudes are environment-averaged values is validated in Section IX by the 
CERN data, validating the view of pion creation processes as open quantum systems interacting 
with environment. In Section X. we introduce a concept of decoherence free subspace that is central 
in determination of the quantum states of the environment, briefly discussed in Section XI. The 
time- irreversibility and CPT violation in pion creation processes is discussed in Section XII. The 
paper closes with conclusions in Section XIII. Appendix A presents the Lutz-Rybicki tables of 
relations for density matrix elements in terms of amplitudes on which the present work is based. 
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II. FINAL STATE DENSITY MATRIX IN tt ~p tt 7r+n AND RECOIL NUCLEON 
POLARIZATION. 

A. Prom S'-matrix to transition matrix T. 

In the laboratory system of the reaction ir~p — ► Tr~ir + n the +z axis has the direction opposite 
to incident pion beam. The +y axis is perpendicular to the scattering plane and has direction of 
Pn- x Ptt-tt+- The angular distribution of produced dipion system is described by the direction of 
7r~ in the dipion center-of-mass system and its solid angle O, = 9, <j). The target nucleon and recoil 
nucleon helicities v and \ are defined in s-channel helicity system. The dipion helicity A will be 



defined in the t-channel helicity system 32|, l4ll . |46| . 1471 ] . 



The pion beam and nucleon target are prepared in a separable state pi(ir p) = pi(ir ) (g) pi(p) 
where Pi(ir~) = \p-kO >< p n 0\ <S> |1, — 1 >< 1, — 1| and 

vv 

IPttO > and \p p v > are pion and proton helicity state vectors, respectively, and \I, 1% >< I, I 3 \ are 

their isospin states; stands for pion helicity. p'p (P) is the target spin density matrix. In the 
following we suppress the momentum and isospin labels in the initial and final helicity states | Ov > 
and \8(j), \ >• The final state vectors form a basis of orthonormal state vectors in the recoil nucleon 
spin space < 9(j),x\6<ft,X >= & xx ' ■ The final state density matrix pf(ir~ir^ 



n) 



p f (9(t>, P) = £>/(00, P ))ll> X >< 0<t>, X\ (2-2) 



xx 



has matrix elements given by evolution equation pj = SpiS + 

= Y, < °^x\S\0u > (pj(P)) w > < Oz/|S+|00,x' > (2.3) 

We will use an abbreviated notation for amplitudes S Xt ov =< x|<S'|0z' > and suppress their 
dependence on energy s, momentum transer t, dipion mass m and the angles 9(f>. We also drop 
the superscripts 1 1 



22' 



Transition T-matrix is defined by Sfi = Ifi + i(2ir) 5 (Pf — Pi)Tfi. With amplitudes S x , 



0;/ 



iT xfiy {2TTf5\P f - Pi) we get 



P f (94>,P) xx ' = Pf{e<t>,P) xx ,{VT){2-KT8\P s - Pi) (2.4) 

where p'^(94>,P) is expressed in terms of transition amplitudes T x ^ u and where we have used the 
conventional approach to deal with a square of (^-functions [3^] with V and T being total volume 
and time confining the interactions to be taken in the limit V, T —* 00. According to the Born rule, 
the probability of n~p(u) — ► tt~ (pi)n + (p2)n(p3, x) is given by 

dP xfiu = \S xM \ 2 [] = \T x ^\ 2 d<S> 3 (P uPl ,p 2 ,p 3 )(VT) (2.5) 

71=1 n 

where the Lorentz invariant phase space d<3?3 = q(m 2 )G(s)drndtdVl with q(m 2 ) and G(s) the pion 
momentum in dipion center-of-mass system and G(s) the energy dependent part of the phase 
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space [52j]. The probability per unit volume, unit time and per target particle is do x $ v = 
dP X) ov/(VTFlux(s)) and the differential cross-section reads 

dcr x ,0u _ q(m 2 )G(s) 2 

dtdmdU Flux(s) 1 x ' 00 ' 1 ' 

Applying formally the same procedure to every bilinear term S > q u 'S x Ql/ of p/(^, P) we can define 
a differential cross-section matrix 

where we have absorbed \/ q(m 2 )G(s) / Flux(s) into transition amplitudes and redefined pf(Q,P). 
It can be written in a matrix form 

Pf = Tp. t T+ (2.8) 

where T is the matrix of transition amplitudes. The transition matrix T is non-unitary and non- 
hermitian but it still evolves pure initial states into pure final states on account of the central 
assumption (1.1) and unitary S'-matrix (1.2). 



B. Conditions for purity of initial and final states. 



The target nucleon spin density matrix has the form 37|, [3 



Pp {P) = \{l + Pa) (2.9) 

where P is the target polarization vector and a = (cr x ,a y ,a z ) are Pauli matrices. The target is 
in a pure state if and only if \P\ 2 = 1; otherwise it is in a mixed state (l3 . 38]. When P = 



the target is maximally mixed state with equal probability of the target nucleon spins "up" or 
"down" relative to the scattering plane, and the spin ensemble is isotropic. When P = (0, ±P, 0), 
the target is transversely polarized with spins "up" for +P or "down" for — P, respectively. For 
P = ±1 the transversely polarized target is in a pure state p P:U = ^(1 + a y ) or p p ^ = |(1 — a y ). 
In modern polarized targets the initial density matrix pi(P) can be varied by external magnetic 
fields to rotate the polarization vector P into any desired direction. 

To discuss the purity of the final state we note a useful result from quantum state tomogra- 
phy [13]. Arbitrary density matrix p of n qubits can be expanded in a form 



r(a Vl (gi a V2 (g) ... (g> a Vn p)a Vl 8) a V2 (8) ... <g> a Vn ( 2 -10) 

where the sum is over the vectors v = (t>i, v%, ■■■,v n ) with entries chosen from the set a J ,j = 0, 1, 2, 3 
of Pauli matrices and a = 1. This result can be generalized to any non-qubit systems The 
traces in (2.10) represent average values of spin correlations. The determination of these averages 
requires repeated measurements forming a large ensemble of events. The expansion (2.10) shows 
that the concept of density matrix in Quantum Theory is inherently a statistical concept. The 
final density matrix pf(9(j), P) is a single qubit density matrix corresponding to spin ^ of the recoil 
nucleon. It can be written in the form (2.10) 

Pf {6<t>,P) = ^(J°(^,P)a° + J>0,P)a) (2.11) 
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where the traces P(9(f),P) = Tr{a 3 pt(9(f),P)),j = 0,1,2,3 represent measurable intensities of 
angular distributions as seen from (2.7). Introducing recoil nucleon polarization vector Q(9(j),P) 
using a relation 

Q(0</>, P)i°(9(f>, P) = T(e<t>, P) (2.12) 

we can write 

p f (9<f>, P) = + Q(9<f>, P)B)I°{94>, P) = p n (Q)I°(9<f>, P) (2.13) 

The normalized final state density matrix pj(9<f>,P) = pf(9(j), P)/I°(9(f), P) is simply the spin 
density matrix of the recoil nucleon p n (Q). It will represent a pure final state if and only if the 
recoil nucleon polarization vector Q = {Q l -,Q 2 -,Q 3 ) satisfies the condition \ Q\ 2 = 1 for all solid 
angles VL = (#,</>) at any given dipion mass m and momentum transfer t lj|, 38]. 

We can see that the vector Q(9<p,P) defined by (2.12) is recoil nucleon polarization from 



the definition of polarization vector [3a . [39J. A spin state of an ensemble of particles with spin 
posseses a vector polarization when in the rest frame of the particle the spin operator s has a 
non-zero expectation value < s >= Tr(sp). In general, a polarization vector Q is defined as 
Q =< s > /(sTr(p)) = Tr(sp) / (sTr(p)) . For the final state ensemble of recoil nucleons s = \a 
and p = pf, and (2.12) holds. 

The polarization vector Q has transverse components Q 2 and Q 1 where Q 2 is perpendicular to 
the scattering plane in the direction of the y axis and Q 1 is transverse to the neutron direction in 
the scattering plane. The longitudinal component Q 3 is along the direction of neutron motion in 
its rest frame. Measurements of Q 2 and Q 1 can be done by rescattering of the recoil neutron on 
spin zero carbon target. Measurements of Q 3 are more difficult since they require spin rotation by 
special magnetic fields to convert the longitudinal polarization to transverse polarization. 

C. Angular expansion of the final state density matrix. 

Using p p = 1(1 + PS) we can write matrix elements of pj in terms of components of target 
polarization 

{pf{H,P))ll = + PMmll; +Py{Py{H))f y! +Pz{Pzm))ll (2-14) 

In (2.14) the subscript u stands for unpolarized target P = 0. Using the decomposition (2.14) of 
Pf(9(j),P) we find a decomposition for the intensities 

P{9<j>,P)=Tr{ai pf {94>,P)) = li{e<l>) + PJi{e<l>) (2.15) 

where the components I J k (9(j)), j = 0, 1, 2, 3 and k = u, x,y, z of the intensities P(9cp, P) are given 
by traces 

ll(9<j>) = Tr XiX ,((^') xx »(p fc (^))j4) (2.16) 

In general, a plane wave helicity state of two particles with helicities pi,p>2 can be expanded in 
terms of angular helicity states 



P1P2 >= E V ^T^ D lv& °' -<j>)\P JX -i > (2.17) 



47T 

J,X 
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where p is the momentum in center-of-mass system and J and A are the two-particle spin and 
helicity, and fi = fii - /i2- For two pions hi = ^ 2 = and i>( (</>, 0, -(/)) = -y/47r/(2 J + 1)Y X J *(0, 4>). 
The final state can be expanded in spherical harmonics 

l^,X>=E y A J *(^)l</A,x> (2.18) 

JA 

where J and A are dipion spin and helicity, respectively. Using (2.18), the angular expansion of 
the transition amplitudes 



• J A 

defines helicity amplitudes of definite dipion spin H XxQu (s,t,m). The partial wave helicity ampli- 
tudes (2.19) are normalized such that the intensity of 7r~7r + production measured on unpolarized 
target and integrated over the angular distribution of the produced pions is given by 

= / dnTr x=x ,(( Pu (e,^,) ^EEKoJ 2 (2-20) 

The target polarization components of matrix elements (2.14) are given by (2.7) 

(P*W))H = \Y, H xM0<t>)(<7k) uu >H* x ,^(94>) (2.21) 



vv 



where k = u,x,y, z and {pv) vv > = S uu ' is a unit matrix. Using (2.19) their angular expansion reads 

(pkM))**, = EE^)jAV y ^^^/*(^) (2-22) 

fx' 

where (^)a^ 'a'x' are ^ e an S u l ar spi n density matrix elements of the final density matrix. In terms 
of partial wave helicity amplitudes 0u these angular matrix elements read 

(R k (s,t,m)) J ^} = IEK,0v(^vv>h{>;>, ; - 0n.{s,t,m))%f}^ (2.23) 

vv 

where we defined normalized angular spin density matrix elements (pk{s, t, m)) 2 ' , 2 , . The compo- 

nent intensities I 3 k (6(p), k = u,x, y, z, j = 0, 1, 2, 3 have angular expansions arising from the traces 
(2.16) 

JA J'A' 

where the unnormalized angular density matrix elements R 3 k are the traces over recoil nucleon 
helicities 

(Ri(s,t,m))({, = Tr x ^(a% x// (R k ) J x lf) x ) (2.25) 
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Expressed in terms of partial wave helicity amplitudes they read 

(^)aa' = (4)av^ = \ E E(^')xx"^A X ",o,(^)-'^Vx,o,' (2-26) 

where we defined normalized angular density matrix elements (pl)\\i ■ The normalization condi- 
tions 

Tr J=f ^{Tr x= ^{p u { S ,t,m)) J ^) = Tr J=y iW ((p£(M,m))$) = 1 (2.27) 

are equivalent to (2.20) for normalized helicity amplitudes. 

Combining in the sum (2.25) the terms with inverted JA and J'A' we can write (2.25) for each 
k = u,y,x,y and j = 0, 1, 2, 3 in the form 

\ E E^aaW Y$* + RLLxY^Y* + R J y J x Yj'Y x J * + R J ^_ X Y^Y^ X ] (2.28) 

JA J'A' 

Using hermiticity of the density matrix 

= (R{) J x J y* (2.29) 

and a relation for spherical harmonics Y^ M (6,(p) = (— 1) M (Y^(9, (j)))* the sum of terms in (2.29) 
takes the form 

[+2Re(R{ J y + (-l) A+A 'i?^'_! A ,) J Re(y/y/*) (2.30) 



-2Im{Rg, + {-l) x+x 'R J j;i x ,)Im{Y x J Y x j ;*)] 



The amplitudes H x 0u describe a two-body process 7r + p — ► " J(m 2 )" + n where J and m 2 are 



the spin and mass of the dipion particle " J(m 2 )". The initial and final states in this process are 
separable. For such states the parity conservation in strong interactions requires that 38, 39. liol] 

^A- x ,o-„ = (-l) X+x+u H J Xxfiu (2.31) 
Parity conservation relations (2.32) imply symmetry relations for spin density matrix elements 

= +(-1) A+V K)- J a'-A' (2-32) 
for (k,j) = (u,0),(y,0),(n,2),(y,2),(x,l),(z,l),(x,3),(z,3) and 

(RDU'i = -(-1) A+A '(i4)-1'-A' (2-33) 

for (x, 0), (z, 0), (x, 2), (z, 2), (u, 1), (y, 1), (u, 3), (y, 3). Using these symmetry relations the compo- 
nents I 3 k (9(j)) of the dipion angular distribution P(6(f),P) measured on polarized target take the 
form 



JA J'A' JA J'A' 



(2.34) 
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for (k,j) = (n,0),(y,0),(n,2),(y,2),(x,l),(z,l),(x,3),(z,3) and 

W) = EEM)av^a(^v*W) = ^EE(W fc )^/m(y/(6i0)y/*(^)) 

JA J'A' JA J'A' 

(2.35) 

for (x, 0), (z, 0), (x, 2), (z, 2), (u, 1), (y, 1), (u, 3), (y, 3). The elements (ImR{)^ in the group (2.34) 
and {ReR 3 k ) J X y in the group (3.35) are not observable as the result of parity conservation in strong 
interactions. However, as we shall see in Section ?? in the case of S- and P-waves, the measured 
elements on polarized targets supplemented by certain phase relations between amplitudes enable 
to calculate the unobservable elements. 



D. Experimental form of angular distributions. 

The expressions (2.34) and (2.35) for angular intensities I J k (0(p) in terms of angular matrix 
elements {RV)\\i can be simplified. Using the hermiticity of the spin density matrix elements 
(Pk)\'\ = (Pk)H'* we can write ( 2 -34) and (2.35) in the form 

J XX' J<J' XX' 

(2.36) 

Iim = ^EE 2( W fc )^,7m(y/(^)y/*(^)) (2.37) 

J<J' XX' 

The number of spin density matrix elements is further reduced by making use of their parity 
symmetry relations (2.32) and (2.33) and relations for spherical harmonics 

Re(Y£ x Y$) = (-l) A+A 'i?e(y/y/*) (2.38) 

Im{Y£ x Y%) = -(-i) A + A 7m(y/y/*) 

We can then simplify the sums 

E = EE+EE = EE&(-) (2.39) 

A, A' A>0 A' A<0 A' A>0 A' 

where £o = 1 an d £a = 2 for A > 0. The two sums in (2.36) can be combined by introducing a 
factor = 1 for J = J' and = 2 for J < J'. The final expressions for (2.34) and (2.35) with 
independent angular density matrix elements then read 

Iim = iiL E EE^<^ Re ^)ii'MYx J (om J *m) (2.40) 

J<J' A>0 A' 

Iim = E EE^W)^^/w)^w)) (2.4i) 

J<J' A>0 A' 



Note that in (2.41) Im{p > k )'l i , l = and Im(YQ \64>)Yq'*{6^>)) = 0. Because of the angular properties 

of Y^(9(j)), the elements (p^)oo = (Pfc)ss' (Pfc)oo and (Pfc)n are n °t independent but appear in two 
independent combinations 

(pDsp - 04)°° + (/4)oo + (p{)pp ^ 04)oo " (pDW (2-42) 
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The most feasible experiments are measurements on unpolarized or polarized targets which 
measure the two-pion angular distribution I°(9<p,P) leaving the recoil nucleon polarization vector 
Q not observed. Such measurements provide information on the reduced final state density matrix 
given by /°(6><£,P) 



P) 



Tr x=x/ ((p f (e<f>,P)^ 



+ P 1° 

I 1 X 1 ' 



+ Pyi$m + p z i° z m 



(2.43) 



It is the intensity I°(9tp, P) which has been measured in CERN measurements of pion creation 
processes on transversely polarized targets and on which this study is based. 



Experimentally, in a given mass region only helicity amplitudes with J < J ma x contribute and 



all sums in (2.40) and (2.41) are finite. The spin observables {pL){{i for J < J max are determined in 



small (m, t) bins from the measured angular distribution P (9(p, P) at any given target polarization 
P using the statistical methods of maximum likelihood function (4ll . 46 . 55j_56]. This process is 
referred to as quantum state tomography in Quantum Information Theory lie 



E. Nucleon helicity and transversity amplitudes with definite f-channel naturality. 



The helicity amplitudes H\ x o„ are combinations of helicity amplitudes with definite t-channel 
naturality r\ = VS where V and 5 are the parity and the signature of Reggeons exchanged in 
7r~ +p —* " J(m 2 )" +n [40j. The natural and unnatural amplitudes 0± and U x+ . corres pon d 
to naturality rj = +1 and rj = — 1, respectively. They are given for A ^ by relations 



u x+,o± 



1 



(hL 0± + (-i) A # 



XttJ 



-x+,o±j 



(2.44) 



(-1) 



X+,0±J 



For A = they are 



J 

0+,0i 



and NL 0± 



(2.45) 



In (2.44) and (2.45) + and - correspond to +5 and 
exchange amplitudes U J 



values of nucleon helicities. The unnatural 
and U£, 0+ exchange it and a\ quantum numbers in the i-channel, 
respectively, while the natural exchange amplitudes N^, q _ and N^, q , both exchange a>2 quantum 
numbers. 



A+,0- 



Amplitude analyses of measurements on polarized tar gets are best performed in terms of 
transversity amplitudes with definite t-channel naturality 



37, ft 



In such measurements 



the spin states of the target nucleon are described by transversity r with r 



+ 



1 



u and 



d corresponding to "up" and "down" orientations of the target spin relative to the 
37 . [3^ . Following Lutz and Rybicki 32], we define mixed helicity-transversity 



scattering plane 

amplitudes with nucleon helicity replaced by nucleon transversity 



u Xt„,0t p 



Xp ,A n 



( 2' 2' 



7T. 
2 



H J 



AA„,0A. 



,7T vr 

v 2' 2' 



7T . 



(2.46) 



where D? 

AT 



,7T 7T 

2' 2' 



7T 



2 )=^ (A - r) i(|) 
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The factor £ = 1 or f = e i7r ( An_A f) for the y axis in the rest frames of target and recoil nucleon in 



the direction or in opposite direction to the normal to the scattering plane, respectively 371.1381. |47I|. 
In our case £ = 1. Then for A = and A / we obtain 

Godflu = ~~ *(-^o+,o+ + iHo+,o-)> Qju.od = +i(Hq +0+ — iHq + q_) (2.47) 
Parity conservation (2.31) requires that 

Qht.Ou = ^Od.Od = G\d,0u = ^Xu,0d = (2.48) 

Omitting the inessential factors ±i in front of parenthesis in (2.47), the unnatural exchange 
transversity amplitudes are given by 

Ul = Ui u = U( Tnfiu = ^(G{ Tnflu + (-l) X G J _ XTnfiu ) = ^(Ul +fl+ + iUt +t0 _) (2.49) 

J 1 1 

= U\,d = U\ Tnt od = —j={G XTn 0d + (-1) G_ XTn 0d ) = -j=(U x+fi+ - iU x+ _) 

where r n = — r p and r n = +r p for A = and A 7^ 0, respectively. The natural exchange transversity 
amplitudes are given by 

N J X = N(, u = N( ufiu = ^=(G J Xufiu - (-l) x G J _ Xufiu ) = ^(^a J +,o+ + i^A+,o-) ( 2 -50) 
N J X = N( !d = N J Xdfid = -^(G( di0d - (-l) x G J _ Xdfid ) = ^(^+,0+ " ^a+,o-) 

iV J = iV J = OoriV ^ = iV ^ = o 
The factor 1 / y/2 results in the normalization 

d2(T f ^(^) = EEi^i 2 + K. p i 2 (2.51) 



JX Tp 



dtdm 

Parity conservation requires 

u j x = +(-irui, u: = +(-iru x (2.52) 

N J = -(-l^Ni, N J X = -(-1) X N J X 

Inverted relations for helicity amplitudes in terms of transversity amplitudes with definite t-channel 
naturality for A > read 



# A J +,o+ = \[{U J X +U J X ) + (N( + N{)\ (2.53) 
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Hi x+ , 0+ = { -^-[(U( + U{) - {N( + N{)\ 



Hl +fl - = -{Wl - U{) + {N( - N{)\ 



Hix+jo- = -A^KU J X - U{) - (N( - N{)] 



General expressions for all spin observables (RV)x\i for several types of helicity and transve rsity 



amplitudes were calculated by Lutz and Rybicki and are given in their unpublished paper Ref. 321 ] . 
Their original Tables la, Id and lb for helicity amplitudes, transversity amplitudes with definite 
t-channel naturality and helicity amplitudes with definite t-channel naturality, respectively, are 
reproduced in the Appendix A. 

F. Trace measurements. 

The traces Tr{Rep' k ) can be measured in experiments on unpolarized or polarized targets with 
or without measurements of recoil nucleon polarization in which the angular distribution of the 
produced pions is not observed. The traces are given by the integrated angular distributions 

J A 

for (k,j) = (u, 0), (y, 0), (u, 2), (y, 2), (x, 1), (z, 1), (x, 3), (z, 3). The integrated angular distributions 

l[ = J ll(9(f>)dn = (2.55) 

vanish for (k,j) = (x, 0), (z, 0), (x, 2), (z, 2), (u, 1), (y, 1), (u, 3), (y, 3) since the integral 

F/(00)y/*(^)dO = 6jj,d xx > (2.56) 

is real |57j. The intensities I®(9(j)) provide information about the polarized target asymetry vector 
T = (T x ,T y ,T z ) defined by 

T fc (^)I°(#) = 4°(^) (2.57) 
The averaged values of its components 

_ fT k (ecb)i° u mdn _ JW)dn 

are T x = T z = and T y = T where T is polarized target asymmetry. With 

d 2 a /" n,„„ ,„ , „ d?a 



[ l9(d</>)dn and T^- = [ I°(6(f>)dn (2.59) 
J dtdm J y 



dtdm 

we get from (2.54) trace conditions 

a = J2(Rep° u )i( = 1 and T = ^(Rep° y )( J x (2.60) 



JX JX 
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When the recoil nucleoli polarization is not measured, the integrated intensity measured on polar- 
ized target 

7V)=7»„ + P 9 7^(l + P S T)-^- (2.61) 
does not depend on components P x and P z of the target polarization. 

When the it~it + or 7r°7r° angular distribution is not observed, we measure the average values 
Q 1 (P) of the components of the recoil nucleon polarization vector 

w [ ' ~ fi°(9ci>)dn ~ fi°(e<f>)<m [ b > 



Using (2.15) we can write 



Q j (P)T° = I j = I 3 U + Pj x + Pyll + P Z I J Z (2.63) 



With averaged intensities (2.55) vanishing due to parity conservation, we have the following result 
for the recoil nucleon polarization four-vector Q J , j = 0, 1, 2, 3 

Q°(P)7°(P) = 7° + P y f y = (Q° u + PyQl)^ (2.64) 
Q 2 (P)t(P) =7 2 u + P y ll = (Q 2 U + P y Q 2 y ) d2a 



dtdm 



}\p)i\p) = pj\ + p z i] = (p x qI + p z q 1 ; d ' n 



dtdm 



q\P)i\P) = pjI + pjI = (p x qI + P z Q 3 z ) d ' 



dtdm 



t2 



For transversely polarized target with P x = P z = only the transverse component Q of the 
averaged recoil nucleon polarization is non-zero. Note that Q (P) = 1. 



G. The structure of the averaged recoil nucleon polarization vector. 

Using the expressions (2.26) and the parity relations (2.31) we can express the traces l\ in 
terms of helicity amplitudes. We find 

t = ^ = E(l^ + ,o + | 2 + (l^ J + ,o-| 2 ) (2.65) 

J A 
J A 
J A 
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J A 



where A and R are spin rotation parameters denned in analogy with irN — * irN scattering 37J, l3a . 



39]. For the remaining traces we get 

T « = E +(-l) A Mtf-A + ,0-tfA J ;,0 + - ^-A+,0+^A J ;,0-) (2-66) 
JA 

A/ = ^2-(- l ) XRe ( H -X+,0+ H X+,0+ + H -X+,0- H X+,0-) 
JA 

JA 

^ = E +(- 1 ) Ai?e (^-A+,0+^A+,0+ + ^-A+,0-^A J ;,0-) 
JA 

In our next step we shall relate the four traces in (2.66) to the four traces in (2.65). After some 
algebraic rearrangements of the terms and sums in (2.66) and using parity relations for the helicity 
amplitudes we find 



dtdm ' 1 dtdm 

For now, the traces a, T, R, A in (2.67) are expressed in terms of helicity amplitudes as shown in 
(2.65). The new terms ojy, tat, Rn, An are expressed in terms of natural exchange helicity am- 
plidudes and are defined next. With A^ = U^, for unnatural and natural exchange amplitudes, 
we define normalized partial wave intensities cr^(JA) and normalized partial wave polarizations 
ta(JX) 

M-tt^lAtf + tfl', r A (J»^\Ai?-ri? (2.08) 
Similarly we define partial wave spin rotation parameters 

Next we need to define their sums for A > 

oa = E °a(JX), t a = t a (JX), (2.70) 

JA>0 JA>0 

A A =Y1 a a(JX), Ra=Y1 R a( jx ) 

JA>0 JA>0 
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When the traces (2.65) and (2.66) are expressed in terms of transversity amplitudes with definite 
naturality using the relations (2.53), we find 

Qu = a u + c;v = a, Q° v = tu + t n = T (2.71) 



2 2 

Q u = ~T + 2t n, Q y = -v + 2crjv 



and 



Q 3 X = A V + A N = A, Q 3 z = Ru + R n = R (2.72) 

Ql = -R + 2R N , Ql = +A-2A N 
where a = 1. The equations (2.72) and (2.73) reveal simple structure of the averaged recoil 

2 

nucleon polarization vector. We make an important observation that the components Q u and 

2 

Q y are fully determined by measurements on polarized target without any measurement of recoil 
nucleon polarization and that these components differ from the traces a and T by natural exchange 
terms on and tn- It is this difference that turned out to be crucial in the initial development of 
experimental tests of unitarity of the ^-matrix. 

III. TEST OF UNITARITY IN irN -> irN AND SIMILAR TWO-BODY PROCESSES. 

The final state density matrix in reaction 7riV — > irN measured on polarized target with polar- 
ization P = {P x , P y ,P z ) has the form 

p f (P) = ±(l°(P) + I(P)a) = ±(l + Q(P)a)l°(P) (3.1) 

where Q(P) is the recoil nucleon polarization. Assuming parity conservation for helicity ampli- 
tudes Ho-A',o-A = — (— 1) A+A Hqx'.ox an d using the definitions (2.26), we obtain the well known 



expressions for components P(P) = Q J (P)I°(P) [30. 1391] 



J°(P) = Q°(P)I°(P) = (1 + TP y )^ (3.2) 



I\P) = Q\P)I\P) = (RP X - AP Z )^ 



I\P) = Q 2 (P)I°(P) = (T + Py 



da 



I\P) = Q 3 (P)/°(P) = (AP X + RP Z )^- 

at 

where Q°(P) = 1. The differential cross section da/dt, polarized target asymmetry T (often called 
polarization P) and spin rotation parameters A and R are defined in terms of helicity or transversity 
amplitudes as follows 

= l-Ho+,0+1 + l^o+,o-| 2 = \H U \ 2 + \Hd\ 2 (3.3) 
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- 2Irn ( H o+,o+HQ + _) - \H U \ 2 - \H d \ 2 
A^ = 2Re(H 0+>0+ H^ _) = 2Im(H d H*) 

= l-^o+, o+| 2 - |#o+,o-| 2 = ^Re(H d H*) 

where the target nucleon transversity amplitudes H T are 

1 1 

H u = -^=(#o+,o+ + ii7o+,o-)) H d = -j=(H 0+t0+ — iiJ +,o-) (3.4) 

The helicity and transversity amplitudes in irN — ► ttN are natural exchange amplitudes. Setting 
unnatural exchange amplitudes in (2.72) and (2.73) equal zero, we recover the relations (3.2). 

The spin observables are not independent but satisfy condition T 2 + A 2 + R 2 = 1 which follows 
from (3.3). Using this condition and assuming that the initial state is pure with \P\ 2 = 1, we 
easily verify that |Q(P)| 2 = 1 and the final state is pure. We also observe from (3.2) that even 
when recoil nucleon polarization is not measured, the measurements of target asymmetry T fully 
determine a single value of the transverse component of recoil nucleon polarization Q 2 . From (3.3) 
we can see that the measurements of T together with da/dt yield a single solution for moduli of 
the transversity amplitudes, as expected from a unitary S-matrix. 



IV. UNITARITY CONDITIONS IN n ~p -> tt 7r+n AND SIMILAR PROCESSES. 

We wish to test unitarity in 7r~p — > 7r~7r + n using exact expressions (2.34) and (2.35) for angular 
components li(0,,P) without truncation of the angular expansion. The final state is pure if and 
only if the recoil nucleon polarization P)\ 2 = 1, or equivalently, if and only if the intensities 

P)\ 2 = (J°(fi, P)) 2 . We recall from (2.15) that P(Sl, P)=I 3 U + P X P X + P y P y + l{P z . We will 
work with two pure initial states with transverse target polarization P = (0, P y , 0) where P y = ±1. 
With these initial pure states the final states will be pure if and only if 

(II ± II? + (4 3 ± II? = (I°u ± I y f ~ ill ± I 2 yf (4-1) 

Using angular expansions (2.34) and (2.35) we can write these unitarity conditions in the form 

(J2 Y,( ImR u ± ImRl) J x J x :im(Y x J Y x j ;*)) 2 + (E E( /mii « ± ImRl){ J x ;irn(Y x J Y x j ;*)) % (4.2) 
JA J'A' JA J'A' 

= (EE( Refl 2 ± ^) A J A j ;fie(y/y/*)) 2 - ± &i?^ A j ;&(y/y/*)) 2 

JA J'A' JA J'A' 

Using the Table l.d of Lutz and Rybicki [32| reproduced in the Appendix A we find for the terms 
with P y = +1 

(ImRl + ImRDU; = -2r ?A 7 ?v i?e(^ A J 7V A / ;* - N(U J X ,*) (4.3) 
(ImRl + ImRl) J x J x : = +2r ]xVx ,Im(U x J N x J !* + N(u{,*) 
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{ReRl + ReR° y )i(' = +2i lxV yRe{U(U J ^ + N J x N J y*) 

(ReR 2 u + ReR 2 y ){{: = -2 w i?e([/ A J ^ A J /* - N J x N J y*) 

where rj\ = 1 for A = and r]\ = l/y/2 for A / 0. The relations for terms with polarization 
P y = — 1 are the same with replacements U X ,N X — ► U X ,N X . In the following we will work with 
the unitarity condition (4.2) for polarization P y = +1 which involves only amplitudes with target 
nucleon transversity "up". The condition for the polarization P y = — 1 will have the same final 
form involving only amplitudes with target nucleon transversity " down" . 

After relabeling, the two terms on l.h.s. of (4.2) have the form 

16R2 = (EE-taWX') - Re{Niui:*))Im(YM'l) 2 (4-4) 
JA J'A' 

JA K/j, 

16/2 = (E E +^xriy(Im(U(N^) + 7m(iV/[//;*))/m(y/y/*)) 2 (4.5) 

J A J' A' 

JA ft> 

Taking squares of the two terms on the r.h.s. of (4.2), the common summation will involve combi- 
nations 

{ReRl + ReR y ){ J y(ReR° u + ReR° y )«*' ~ + ReR 2 y ) J x {, (ReR 2 u + fleflj)*?' = (4.6) 

4r ?A r ?A ,r ?/lV (2i?e(C//C/ A j ;*)i?e(iV / fiV^'*) + 2Re(N(N(!*)Re(U^U^*)) 

After relabeling the second term on r.h.s. of (4.6), the r.h.s. of unitarity condition (4.2) factorizes 
into a product 

16UN = 1S(Y,Y, Re(U-M*)Re{YM'l) (E E R<N* N«'*)Re{Y£ Y* *)) (4.7) 

JA J'A' Kji K'jx' 

Using parity relations (5.52) for transversity amplitudes and the fact that Nq = 0, we can rewrite 
the expression for the sums R, I, U and N in the form 

R = -2 E E ^V,Re{U(N^)Re{Y x J )Im{Y^) (4.8) 

J,\>0K,fi>0 

' = -2 E E ^nx^Im{UiN^)IU{Y^)Im{Y^) 

J,A>0 A",/i>0 



20 



j,a>o j',y>o 



^ = 4 E E r,Mte(NiNX'*)Im{YF)Im(Y«') 

K,n>0 K',n'>0 

where £\ = 1 for A = and ^ = 2 for A > 0. The unitarity condition R 2 + I 2 = UN then reads 
E E E E C xx ,^,(Re(U J x N«*)Re(U J x !N«'*)+ (4.9) 

J,A>0 J',X'>0 K,n>0 K',fj,'>0 

Im(U(N^)Im{u{,N^*))Y^^ K '(n) 

= E E E E c^(Muiu^)MN^N^))Y^ K \a) 

J,A>0 J',A'>0 K",/i>0 K',n'>0 

where 

n J A',y («) = i? e (^)i?e(y/)7m(y^)/m(yj) 

The condition (4.9) must hold true for all solid angles fi. This is possible if and only if 

Re(U(Nf*)Re(U J x ;N«'*) + Im(U J x N^)Im(U J x ! N^'*) = Re(U(Uy*)Re(N^* N^'*) (4.11) 

is true for all J, J', K, K', A, A', fi, fi'. With target polarization P y = —1 we obtain from (4.2) 

the same unitarity condition (4.11) for transversity amplitudes U x and with target nucleon 
transversity " down" . Next we introduce relative phases 

«(a J ' = «^a J ) " P«f = *«) - (4-12) 

7a J ^ = <^a J ) " *(^), 7a J 'T = ) " *(JVf) 

and similar relative phases a X y, , j x Jf and 7^'^' for transversity amplitudes U x and . 
With ' 

JK J'K' _ n J J' oKK' —JK -J'K' _ —J J' -nKK' , . 

7a m -7av -«aa'-/V . 7a m -7av - «aa' - V I 4 - 13 ) 
the two unitarity conditions (4.11) now have a final form 

|C//||C/ A J /|Kl|iVf |cos(a^' -d K J) = \U(\\U(;\\N^\\N^'\cos(ai{;)cos(^f) (4.14) 

l^lFi'll^ll^'lcos^;) --ft* ) = |Z7^ 1 1 17^' 1 1 iV^ 1 1 iV^' I cos («^' ) cos (^^' ) 

The unitarity conditions (4.14) will be satisfied in two cases. In the case that all moduli are 
non-zero, the conditions (4.14) require 

sin(a^) sm(^f) = 0, sin(aft') sin(^f ) = (4.15) 
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The conditions (4.14) will also be satisfied when only one unnatural or only one natural exchange 
amplitude is non-zero. In such 

\U(\ = \U J x: \SjjJxx u , \U J x \ = \U J x d d \5jj d 5 XXd or (4.16) 



Kl = \N^\S KK J^, \N„ | = \N/ d \5 KKd 5^ d 

Unitarity of the 5-matrix thus requires that in n~p — > ir~ir + n and other similar processes either 
the condition (4.15) (requirement A) or the condition (4.16) (requirement B) be satisfied for all 
J, J' , K, K' , A, A', /j,, // at all energies s, momentum transfers t and dimeson masses m. 



V. TESTS OF UNITARITY CONDITIONS IN MEASUREMENTS OF tt p -> tt 7r+n ON 
POLARIZED TARGET. 

To discuss the experimental results we first introduce the spectroscopic notation for the ampli- 
tudes 

S T = £/° T , P T ° = 0* T> P? = Ul T , P T N = Nl T (5.1) 



tt 2 n u 



uL,d? 



Nl T ,Df = Ul T ,D™ 



where r is the target transversity r = u, d. The CERN measurements and amplitude analysis of 
ir~p — * Tr^Tr + n on polarized target at 17.2 GeV/c fall into 3 kinematic regions. 

In the first region dipion masses are in the range of 580-900 MeV at low \t\ (0.005 < \t\ < 0.2 
(GeV/c) 2 ) and in the \t\ range of 0.0 < |i| < 1.0 (GeV/c) 2 at p°(770) 



41, 51 



mass. In both these 
[52l . liil all agree on 



subregions S- and P-wave dominate. Amplitude analyses of this data 
non-zero moduli for all amplitudes and non-zero relative phases between all unnatural exchange 
amplitudes. Since there is only one natural this data do not violate unitarity conditions (4.15) and 
(4.16). Explicit calculation shows that |Q(f2, P)\ 2 = 1 for pure states P = (0, ±1, 0) for all solutions 
of amplitudes. Similar results were obtained in the amplitude analyses of CERN measurements 

of 7T 



n 



7T 1 7T p 



46j, |4g, |51|, |52j] at 5.98 and 11.85 GeV/c, and K~ 



n 



P 



43, |5flJ at 5.98 



GeV/c for dimeson masses below 1000 MeV. 



In the second region the dipion masses extend from 580 - 1780 MeV at low momentum 
transfers 0.01 < \t\ < 0.20 (GeV/c) 2 . Two amplitude analyses of the data with different 



mass binnings [42J, |43j] found non-zero moduli of all amplitudes with helicities less than 2, and 



reported non-zero relative phases between unnatural echange amplitudes. There are two pairs of 
natural exchange amplitudes P^ and , r = u,d whose relative phases &(P N ) — $(D Ar ) are 
not reported. Unitarity condition (4.15) requires that the phases (3{l = $(P*) - &(D*) and 



$(Pi^) — $(D^) vanish. The data violate the unitarity condition (5.16) but in the absence 



of data on &(P ) — <fr(D ) the unitarity condition (4.15) cannot be tested in this kinematic region. 



In the third region the dipion masses are in the range 580-1480 MeV at large momentum 
transfers 0.2 < \t\ < 1.0 (GeV/c) 2 . The analysis of Rybicki and Sakrejda of this data 44] found 
non-zero moduli of all amplitudes including helicity 2 amplitudes D 2U and D 2N , a finding which 
violates unitarity condition (4.16). In addition to non-zero relative phases between unnatural 
amplitudes they report a non-vanishing phases $(P^) — ^(D^) and $(Pj^) — §(Dj ), a finding 
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FIG. 1: Relative phases $(P^) - $(-Df ) and - $(£>?). Data from Rybicki and Sakrejda [44| 
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FIG. 2: Relative phases $(P r °) - and - <S>(Df). Data from Rybicki and Sakrejda 
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that violates the unitarity condition (4.15). On the basis of this evidence we conclude that pure 
states can evolve into mixed states in n~p — ► ir~Tr + n and that unitarity is violated. 

Figures 1 and 2 show selected relative phases from the analysis of Rybicki and Sakre- 
jda. The phases $(-0 - $(P>^) and - $(D$) are shown in Figure 1. Phases 

§{D* V ) - $(D°),$(.P°) - $(P>°) and $(P r °) - ^(Df) shown in Figures 1 and 2 illustrate the 
relative phases between unnutural exchage amplitudes. They will be used in Section VIII. to 
confirm a crucial prediction arising from the principal consequence of violation of unitarity, 
namely that pion creation processes behave as open quantum systems interacting with a quantum 
environment. 

We close this discussion with a comment on BNL measurement of it~p — ► 7r°7r°n at 18.3 GeV/c 
on unpolarized target [45| . In this high statistics experiment the dipion mass is in the range 350- 
2200 MeV for several binnings of momentum transfers at low and high \t\. The analysis is made 
on assumption that amplitudes with dipion helicity 2 do not contribute in D-wave and amplitudes 
with nonzero dipion helicity do not contribute in G-wave. When this assumption is removed the 
analysis yields non-zero amplitudes in all waves, albeit with ambiguities that cannot be resolved 
by the measurements on unpolarized target. Measurements of this process on polarized target [33| 
may resolve these ambiguities and thus provide independent evidence for violation of unitarity in 
pion creation processes. 



VI. UNITARITY, SUPERSCATTERING OPERATORS AND THE METRIC. 

With the spinor representations (2.9) and (2.13) of the target and recoil nucleon spin density 
matrices pi{P) = P P (P) and Pf(P) = p n (Q(P))l (P) we can associate four-vectors P M = (1,P) 
and Q a I° = (1,Q)I°, respectively, representing the polarization four-vectors in the rest frame of 
the particles. The determinants of the density matrices are the Lorentz invariant norm of the 
four-vectors 

det(p p (P)) = 1 - |P| 2 = v^P^P" = \\P\\ 2 (6.1) 



det(p n (Q))(I ) 2 = (1 - |Q| 2 )(/ ) 2 = r) a/3 Q a QP(I ) 2 = ||Q|| 2 (/ ) 2 

where r)uv = diag(+l, —1, —1, —1) is metric in Minkovski space-time. The target and recoil nucleon 
polarization four-vectors are connected by the equations (2.15) 

Q a {VL, P)J°(0, P) = I^{VL)P^ (6.2) 

In matrix form we can write QI° = IP so that the matrix I has the meaning of Hawking's 
super scattering operator connecting initial and final states 0]. The norm ||Q|| 2 (^ ) 2 now takes the 
form 

HQIIV ) 2 = Va/3^^P^P U = G^(n)P»P" (6.3) 

Unitary 5-matrix requires that any pure state ||P|| 2 = evolves into a pure state IIQH 2 = 0. A 
sufficient condition for this to happen is for the correlation matrix 

G^(O) = Z(Q) Vflv (6.4) 
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for then ||Q|| 2 (/ ) 2 = Z(n)\\P\\ 2 
the Minkovski metric in (6.4) we turn to ttN 
(3.2) with the matrix If* given by 



0. To find out what could be the matrix Iff that leaves invariant 



ttN. The transformation (6.2) now has the form 



/ 1 T \ 

do R -A 

~dt T 1 

\ A R j 



(6.5) 



The matrix (6.5) is orthogonal but not unitary: 
7-2 T2 



while V J«I, 



/9 



for a ^ /3, 

E M « = T^Ipl = C 1 + ^){da/dtf and £^ = £ M J*i* = (1 - T 2 )(do/dtf with 
determinant equal to (1 — T 2 ) 2 (do / dt) A . We easily verify that the matrix (6.5) preserves the 



Minkovski metric rj a/3 I°I^ = (1 - T 2 )(do / dtfr]^. 



If we impose the structure of matrix (6.5) on the matix i^(f2) ^ n (6-2) 



we obtain constraints on 





its elements which imply constraints on the amplitudes. Prom the requirements that J* = I* 
and J-j = II = we obtain, using the Table l.d of Lutz-Rybicki in Appendix A, conditions 



ReiUlN^- 



■NiU$*- 



-uW'r+NW'*) 



Re(U(N(!*- 



-nzuz;*+u j x n j x*-n j x u j x*) 



(6.6) 



Im(U J x N J x !* + N(U J X ;* + U J X N J X :* + ivX*) = Im(U^!* + JV/tT^ - C/X'* " ^5*) = 

From these conditions we find 

Re(N(ui'*) = Im{N(ui'*) = 0, Re(N J x U J '*) = Im(N J x U J '*) = (6.7) 

These conditions require that either all unnatural exchange ampliyudes with zero helicity or all 
natural exchange amplitudes vanish, in a clear contradiction with known data discussed in the 
previous Section. 

Based on the violation of unitarity conditions for the special transverse initial polarization 
P y = dbl, we conjecture that in pion creation processes the superscattering operatotors in (6.2) 
do not preserve Minkovski metric and thus evolve any initial pure state into a mixed final state. 
Formally, the particle infraction results in a non-unitary transformation of the four-vector 
into the four- vector Q a which transforms Minkovski metric T}^ into metric G^ V (Q) in a curved 
space-time associated with the interaction process. 



VII. NON-UNITARY EVOLUTION OF OPEN QUANTUM SYSTEMS. 

In this Section we briefly review how evolution from pure states into mixed states occurs in 
the theory of open quantum systems interacting with a quantum environment. These concepts 
are essential for the physical interpretation of the violations of unitarity conditions and the 
reinterpretation of the spin formalism developed in the Section II. to form a new picture of 
pion creation processes as open quantum systems interacting with a quantum environment. The 
discussion follows largely the book by Nielsen and Chuang [131 ] . 

Associated with any Hilbert space H of state vectors with basis states \m > is a Hilbert- 
Schmidt space 13(H) of density matrix operators with basis operators \m >< n\ and inner product 
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(A, B) = Tr(A+B) 0, EIQ. Density matrix operator p represents the quantum state of an 



quantum system. It is a positive definite hermitian matrix that satisfies a condition 12, 13, 3j 

Tr(p 2 ) < {Tr{p)f (7.1) 

A state about which we can have a complete knowledge is said to be a pure state. A state of which 
we can have only an incomplete knowledge is called mixed state. The equality Tr(p 2 ) = (Tr(p)) 2 
is satisfied if and only if the state is pure. A quantum state is represented also by a state vector 
\fy > in a Hilbert space H if and only if the state p is pure. In that case p = |^ >< Vl/|. The 
Hilbert-Schmidt space B{H) includes mixed states to which correspond no vector states in H. 

Unitary S"-matrix connects initial vector states \ipi n > in Hi n to final vector states \ipout > 
in H out and the corresponding pure states pi n and p ou t in B(H; Lri ) and B(H out ), respectively. It 
also connects mixed initial states pi n with mixed final states p ou t since such states are incoherent 
sum of corresponding pure states. The question arises whether a linear mapping from B(Hi n ) 
to B(H out ) exists that preserves superposition principle but is more general than the 5-matrix, 
and what physical reality it would represent. Such mappings in fact exist and describe generally 
non- unitary evolution P nn.ti S) = £{pi n {S)) of an open quantum system or process S interacting 
with an environment E 11, 3, 

The co-evolution of an open quantum system S and a quantum environment E is assumed to be 
described by a unitary operator U. The initial state pi(S,E) of the combined system is prepared 
in a separable state pi(S,E) = pi(S) Pi(E). A unitary quantum operation £ describing the 
interaction of S and E evolves Pi(S, E) into a final state [ill , lil . 14 ] 



p f (S,E) = £{ Pi {S,E)) = U Pi (S,E)U+ (7.2) 

Pf(S,E) is not a separable but an entangled state of S and E. The observer can perform mea- 
surements on the system S but cannot perform direct measurements on the environment E. After 
the transformation U, the system S no longer interacts with environment E. The quantum state 
of system S is then fully described by reduced density matrix 

p f (S)=Tr E (pf(S,E))=£( Pi (S)) (7.3) 

in a sense that we can calculate average values < Q >= Tr(Qpf(S)) of any observable Q. The 
trace in (7.3) is over the interacting degrees of freedom \eg > of the environment which form an 
orthonormal basis in the finite dimensional Hilbert space H(E) of the environment E. In general, 
the reduced state Pf(S) is no longer related to the initial state Pi(S) by a unitary transformation 
Pf = SpiS + . Instead it is given by the non-unitary Kraus representation 11, 13|, 14] that generalizes 
the evolution equation pj = Spi,S + 

Pf (S) = J2SePi(S)S+ (7.4) 

t 

where St =< ei\U\eo > are operators acting on the state space of the system S and where we 
assumed that Pi(E) = |eo >< eo| is a pure state. Physically Kraus representation corresponds to 
a measurement of the environment E by the system S just after the unitary co-evolution U. The 
operator elements Si must satisfy the completness relation 



St = I (7.5) 
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for trace preserving quantum operations £(pi(S, E)). When the initial state of the environment is 
a mixed state pi (E) = ]T p mn \e m >< e n \, the Kraus representation takes the form [lit ] 

m,n 

pf( s ^ = yzyzp mnSe - m P^ s ) s ni ( 7 - 6 ) 

i m,n 

where the operators St m =< eg\U\e m > satisfy a completness relation similar to (7.5) and 

Tr( Pi (E)) = Y^Pmm = 1 (7.7) 

m 

When certain degrees of freedom S of the system S are not measured, the density matrix is further 
reduced by taking a trace over S . In certain cases the resulting reduced density matrix Pf(S ) has 
a part in the unitary form Pf = SpiS + independent of the interacting degrees of freedom of the 
environment. This part is referred to as a decoherence free subspace of Hilbert-Schmidt space B(H). 

The Hilbert space H{E) of the environment E is formed by the interacting degrees of freedom 
\e£ > involved in its interaction with the system or process 5. The Hilbert space H(E) has a finite 
dimension dimH(E) < dim Hi (S) dim Hf(S) [ljj. Interactions of S with E in which there is an 
exchange of energy and/or momentum are called dissipative interactions. In contrast, there is no 
exchange of energy or momentun in non-dissipative interactions of 5 with E . Non-dissipative 
interactions are also referred to as dephasing interactions since in general they effect only a change 
of phases. 



VIII. PION CREATION PROCESSES AS OPEN QUANTUM SYSTEMS INTERACTING 
WITH A QUANTUM ENVIRONMENT. 

Assuming unitary S'-matrix and that the pion creation process n~p — > Tr~ir + n behaves as an 
isolated quantum system S, we used evolution equation 

Pf (S) = S Pi (S)S + (8.1) 

to develop in Section II. a spin formalism to express Pf(S) = pf (O, P) in terms of measurable spin 
density matrix elements {RVjiy which are bilinear combinations of transversity amplitudes with 
definite f-channel naturality U^ T and . In developing this formalism we used the form (8.1) 
but we made no explicit use of the unitarity of the 5-matrix or of the assuption that the pion 
creation process is an isolated event. 

Unitary 5-matrix evolves isolated system in any pure state into an isolated system in a pure 
state. In Section IV. we imposed this condition on pf(Tl,P) for specific initial pure states with 
nucleon polarization P = (0, ±1,0) and obtained unitarity conditions (4.15) and (4.16) that are 
violated by the CERN data on polarized target. We concluded that pure inital states p% evolve 
into mixed final states pf in pion creation processes. This means the evolution in these processes 
is not unitary. In quantum theory such non-unitary evolution occurs in open quantum systems 
interacting with an environment. We are thus led to the conclusion that pion creation processes 
are behaving as open quantum systems S interacting with a quantum environment £ in a unitary 
co-evolution described by a co-evolution operator U. The observed final state Pf(S) is then reduced 
density matrix given by Kraus representation (7.6) 

Pf(S) = Tr E (pf(S,E)) = ^2^2pmnSlrnPi{S)S+ e = PmnPf (^71, n€) (8.2) 

I m,n £ m.n 
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To summarize, on one hand we have the spin formalism developed in Section II. from the 
unitary evolution equation (8.1) that has been used in amplitude analyses of the data. On the 
other hand, the data analyses lead to the conclusion that pf is given by a non- unitary Kraus 
representation (8.2). Self-consistency requires that Kraus representation leaves invariant the 
experimental form of angular intensities given by (2.40) and (2.41) and the form of equations for 
density matrix elements in terms of amplitudes given by the Lutz-Rybicki tables. In this Section 
we show that such invariance holds provided that the unitary co-evolution conserves -P-parity and 
quantum numbers of the environment. 

Each term in Kraus representation (8.2) of the mixed state Pf(S) has the formal form of the 
evolution equation (8.1). Thus we can apply the spin formalism of the Section II. to each term 
Pf(£m,n£) = S£ m pi(S)S^ separately with separate co-evolution helicity amplitudes 

H( xfiu (£m) =< JX, x \S tm \0u >=< JX, X \ < e e \U\e m > \0u > (8.3) 

and the corresponding co-evolution transversity amplitudes U^ T (£m) and N^ T {£m). Each term in 
(8.2) has the form (2.11) 

p f (9<t>,P,£m,n£) = \{l\e<p,P) lm)nl a Q + I(9<t>, P) tm>n t3) (8.4) 
with a decomposition (2.15) for the intensities 

P(9<p,P) emM = Tr(aip f (6<f>,P,£m,n£)) (8.5) 

= Ii{04>)tm,nl + PxI J x {04>)tm,nt + PyP y {H) M + P Z P z {H)tm,nt 
The polarization components of intensities I 3 k {04>)tm,nl have angular expansion similar to (2.24) 
4W)tm, n e = Y,Y.W m ^))xw Y x{H)Yi:*{e<j>) (8.6) 

JX J'X' 

where the unnormalized angular density matrix elements R 3 k (£m,n£) are expressed in terms of 
co-evolution helicity amplitudes by relations similar to (2.26) 

(R{(£m, n£)) J x {; = ± £ E(^)xx»4">N(^)^<o^K) (8-7) 
x,x" w 

As the result of linearity of the Kraus representation, the observed final state density matrix still 
has the form (2.11) 

Pf (9<f>,P) = ^{l (9<P,P)a° + I(9<P,P)a) (8.8) 

where 

P k (9ct>, P) = J2 E*W( 7 W) W (8.9) 

£ mn 

The observed intensities P k (9(j),P) are co-evolution intensities li(9((), P)im,n£ averaged over 
the initial state of the environment and summed over its final states. There are four in- 
teracting degrees of freedom of the environment \e£ >,£ = 1,4 allowed by the condition 
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dim H(E) < dim Hi(S) dim H f (S) = {2s p + l)(2s„ + 1) =4. 

To bring If, (6(f), P) to the form (2.40) and (2.41), we need to bring to this form the co-evolution 
intensities I J k (0(j), P)e m ,ne- To this end we follow the proceedure of Section II. C. First we arrange 
the sum (8.6) in the form (2.28). To bring this form to the form (2.30) we need a hermiticity 
relation 

(Ri(£m, nt))%i = (R{(£m, (8.10) 
The actual calculation using (8.7) gives 

(Ri(em,ne)) J y J x = (Ri(n£,£m)) J x {,* (8.11) 

The hermiticity condition (8.8) can be satisfied only when the amplitudes have a diagonal form 

HlxflAZm) = H( xfiv {££)5 lm = H( xflv (l)5 im (8.12) 

Physically this means that the quantum numbers of the quantum states >,£ = 1,4 must be 
conserved in the unitary co-evolution process U. 

To proceed next, the co-evolution density matrix elements R k (£,£))y^, must satisfy parity re- 
lations (2.32) and (2.33). This requires that P-parity is conserved in the co-evolution process so 
that the co-evolution amplitudes satisfy parity relations similar to (2.31) 

Hi x _ xfi _ v {£) = (-l) x+ ^H{ xfiu (£) (8.13) 

With diagonal form of amplitudes (8.12) the components of intensities I J k (6<p)£m,n£ have a diagonal 
form 

I{ (9<f>)tm,rU = A (0<f>)tAmSnt (8-14) 

The diagonal components I 3 k (d(j>)(.^ now have the desired form (2.40) 

limu = E EE^'^(^i(^^))AA'^(y/(^)y^(^)) (8.i5) 

J<J' A>0 A' 

for (k,j) = (n,0),(y,0),(n,2),(y,2),(x,l),(z,l),(x,3),(z,3) and the form (2.41) 

imke = E EE^'^( 7 ^i^^))^'^(^(^)^'*(^)) (8.i6) 

J<J' A>0 A' 

for (x,0), (z, 0), (a;, 2), (z, 2), (u, 1), (y, 1), (u, 3), (y, 3). The measured intensities 11(00) will thus 
have the same angular expansions (2.40) 

W) = E EE^N)M(*'W) (s- 17 ) 

J<J' A>0 A' 

for (k,j) = (n,0),(y,0),(n,2),(y,2),(x,l),(z,l),(x,3),(z,3) and (2.41) 

W) = E EE^'^( /m ^)AA'^(n J (^)^'*(^)) (8.i8) 

J<J'A>0 A' 
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(Rfyxx* in (8.17) and (8.18) are environment-averaged co-evolution density matrix elements 



for {x, 0), (z, 0), (x, 2), (z, 2), (u, 1), (y, 1), (u, 3), (y, 3). The measured density matrix elements 

(4)aa'=Ep«(^(^0)av (8- 19 ) 

The co-evolution density matrix elements (R 3 k (t, are sums of bilinear terms of co-evolution 

amplitudes A J ^ T (£)B^'*(£) given in the Lutz-Rybicki tables in Appendix A. The co-evolution 
transversity amplitudes 

A J ^ T {t) =< J\r),T n \ < ei\U\ei > |0r > (8.20) 



= -j=(< JA,r n | < ei\U\ei > |0r > -??(-l) A < J - A, r n | < e/|f7|ej > |0r >) 

describe the interaction with the environment. Recall that recoil transversity r n is uniquely deter- 
mined by target transversity r and naturality rj = ±1 as the result of parity conservation. In the 
notation of Section II 

U{ T {1) = A J £\£), N{ T {i) = A J £\£) (8.21) 

The measured density matrix elements (R J k )\\i given by (8.19) will keep the form given by Lutz- 
Rybicki table l.d with moduli and bilinear terms of transversity amplitudes A^ v r now redefined 
as environment-averaged moduli and bilinear terms of co-evolution amplitudes, respectively. The 
equations for the measured moduli of transversity amplitudes then read 

|4*;| 2 = Ep«I<?W! 2 > r = u,d (8.22) 
t=\ 

The measured bilinear terms have the form 

Be{A*»X>$*) - \ A \t\\ B \'^*\ c°s(*(4X'?*)) = jlPuRe{A^{i)B^;{l)) (8.23) 



The observed final state density matrix pf thus has the form identical to the form developed in 
Section II. from the evolution equation (8.1). However, it is now the averaged density matrix 
elements (8.19) and the averaged moduli (8.22) and averaged correlations (8.23) that are measured 
in experiments on polarized targets and on which we may no longer impose unitarity conditions. 



IX. EXPERIMENTAL TEST OF KRAUS REPRESETATION FOR OBSERVED FINAL 
STATE IN n-p -> ir-ir+n. 

We have arrived at a crucial prediction arising from the principal consequence of violation of 
unitarity, namely that pion creation process behaves as an open quantum system interacting with 
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FIG. 3: Phase gap A T and cosine gap T T for aplitudes P2, D° and D^F at large momentum transfers 
0.2 < \t\ < 1.0 (GeV/c) 2 . Data from Rybicki and Sakrejda jSj- 

an environment with final states given by Kraus representation. The prediction that the measured 
bilinear terms are environment-averaged bilinear terms of co-evolution amplitudes is testable using 
already existing data on n~p — > 7r~ir + n at large dipion masses. 

The CERN measurements on transversely polarized targets measure correlations of the 
type Re(AB*) = \A\\B\ cos(&(AB*)) where cos(<3?(/4P*)) is the correlation factor. If unitarity 
holds, all these correlations define relative phases <fr(AB*) = &(A) — &(B) between any pair 
of complex amplitudes A and B since there are no averaged bilinear terms. In contrast, the 
measured averaged bilinear terms (8.22) and (8.23) cannot, in general, be represented as real 
parts of bilinear products AB* of two complex amplitudes A and B. We can see that from the 
fact that the averaged bilinear terms (8.23) do not satisfy, in general, the necessary condition 
(Re(AB*)) 2 + (Im(AB*)) 2 = |^| 2 |_B| 2 . This difference allows us to test the validity of the Kraus 
representation. 

We will focus on three amplitudes P°, D® and D 211 . If unitarity holds, their relative phases 
must satisfy phase conditions 

- ($(P T °) - $(£>?)) + ($(P r °) - <5>(D 2 T D )) + (${Dl D ) - = (9.1) 

or, equivalently, cosine conditions for the measured cosines 

cos($(P r °L>°*)) 2 + cos($(P r °D 2l/ *)) 2 + cos($(L> 2f/ P°*)) 2 (9.2) 

-2cos(${P?D T *))cos{$(P?D 2 T u *))cos(${D 2U D° T *)) = 1 
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If unitarity is violated, the measurement of averaged bilinear terms will lead to a phase gap 

- $(P r °Z)°*) + $(P?D 2U *) + §{DfDl*) = A T + (9.3) 

and to a cosine gap 

cos($(P T °I^*)) 2 + cos{^(P^D 2 T u *)) 2 + cos{$(D 2U D° T *)) 2 (9.4) 

-2cos{$(P?D T *))cos(${P?D 2U *))cos(^D 2 T u D° T *)) = T T ^l 
The phases &(P®D®*), &(P®D 2U *) and <$>(D 2U D®*) were measured at large momentum trans- 



fers [44| and are shown in Figures 1 and 2. The phase gaps A T and cosine gaps T T are shown 
in Figure 3. The cosine gaps were calculated only for the mean values of the phases as it is 
very difficult to estimate errors on r r due to nonlinearity in cosines. The phase gaps show large 
deviations from unitarity value 0. The cosine gaps show large deviations from unitarity value 1. 



Both tests are consistent with the observation that unitarity conditions (4.15) an (4.16) are 
violated by the CERN data. Together these facts establish the validity of Kraus representation of 
the observed final states in pion creation processes. They lead us to a new view of these processes 
as open quantum systems interacting with a quantum environment. 



X. DECOHERENCE FREE SUBSPACE. 



In this Section we introduce important concepts of decoherence subspace and decoherence free 
subspace of the Hilbert-Schmidt space of density matrix elements. The concept of decoherence 
free subspace will be cetral to the determination of quantum states of the environment in a sequel 



to this work 60]. 



Experiments on transversely polarized targets with polarization P = (0, P y , 0) measure angular 
distributions I°(Q,P y ). In terms of (2.15) and (2.40) they have a form 

i (n,p y ) = i° u (n) + p y i° y (n)= (10.1) 



E EE^'^((^<)aA' +P y (ReR°) J x (,)Re(Y x J (n)Yt*(n)) 

J<J> A>0 A' 

Maximum likelihood fits to I°(Q,P y ) with J < Jmaxim) determine density matrix elements 
{ReR^) xx , and (ReR° y ) xx ,, or equivalently, the density matrix elements for pure initial states 
P y = P T = ±1, r = u,d 

(ReR (P T )) J x J x , = l -{{ReRl) J x {, + P T (ReR y ) J x J x ,) (10.2) 

It follows from Lutz-Rybicki Table l.d in Appendix A that density matrix elements (ReR° (P T )) XX , 
are related only to amplitudes with the same transversity r = u or r = d. The bilinear terms of 
amplitudes P®, D® and D 2U discussed in the previous Section are in fact density matrix elements 

Re(P?D° T *) = (ReR (P T ))l 2 = \ 2>*((ltefl5(/,/))g + P T {ReR% (10.3) 
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Re(P?D 2 T u *) = (ReR°(P T ))H = -J2Pu((ReR° u (£,£))ll + Pr(ReR y (£,£))H) 

e=i 

i 4 

Re(D 2U D° T *) = (ReR°(P T ))f = - J>«((i?e<(^))i + P T (ReR° y (£, £))%) 

where we made use of relations (8.21) to express the measured elements in terms of the co-evolution 
elements ReR° u (£,£) and ReR° y (£,£). Comp aring with (8.20), we can express the co-evolution 
elements in terms of co-evolution amplitudes 

(ReR° u (£,e))H + P T (ReR° y (£,£))l 2 = Re(P?(£)D° T *(£)) (10.4) 
(ReR° u (£,£))H + P T (ReR° y (£,£))H = Re(P?(l)D 2 T u *(£)) 

(ReR° u (£,£))f + P T (ReR y (£,£))H = Re(D 2 T u (£) D° T * (£)) 

When the bilinear terms (ReA x (£)By*(£)) for a pair of co-evolution amplitudes Af(£) and 
By (£) are different for different interacting degrees of freedom of the environment £, so will 
be the corresponding co-evolution density matrix elements (ReR^(£,£)) X y + P T (ReR y (£,£)) X y . 
As the result, the averaged bilinear term Re(A XT B x '*) does not represent a bilinear product 
of two complex functions and the correlation angle &(A{ Bj'*) does not represent a relative 
phase. We say the corresponding density matrix elements belong to decoherence subspace of the 
Hilbert-Schmidt space of density matrices. 

Consider the case of co-evolution amplitudes -P^(^), D®(£) and D 21, '(£)■ They are complex 
functions and their relative phases must satisfy the phase condition (9.1) and cosine condition 
(9.2). These conditions are violated by the the measured correlation phases from which we 
conclude that some of the measured density matrix elements in (10.3) are environment dependent 
and thus belong to the decoherence subspace. 

Decoherence free subspace is a subspace of density matrix elements that do not depend on the 
interaction with the environment. A density matrix element 

4 

(ReR (P T )) J x J y = Re{A J XT By* T ) = J2puRe(Ai :T (£)B J y* T (£)) (10.5) 

will be environment independent only when the bilinear terms 

Re{Ai T {£)B J y* T {£)) = \A j Kt {£)\\ B J y* T {£)\ cos{<S>{A j Kt {£)) - *(< T W)) (10-6) 

of its co-evolution amplitudes are all equal and thus do not depend on the environment degrees of 
freedom £. In such a case (10.5) takes the form 

(ReR°(P T ))( J y = Re(A J KT B J y* T ) = Re(A((£)B J y*(£)) (10.7) 

4 

where we used pee = 1- The relation (10.7) holds for all £. Note that the moduli and relative 
phases in (10.6) will, in general, still depend on £. The measured bilinear term Re(A x T By* T ) now 
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FIG. 4: Phase gap A r and cosine gap T T for aplitudes P°, D® and S T at small momentum transfers 
0.01 < |i| < 0.2 (GeV/c) 2 . Data from Becker et al. H|. 



has the form of a bilinear product of two complex functions with a well defined relative phase, 
albeit not necessarily a unique one. 

Of central importance is the decoherence free subspace formed by the S- and P-wave subspace 
of reduced density matrix p°f(Q,P) [60.]. In this case the trio of amplitudes S T ,P®,P^ satisfies 
the phase and cosine conditions (9.1) and (9.2 ) which renders the system of equations for moduli 
and phases analytically solvable. In Ref. [60J we show that the quantum numbers labeling the 
solutions can be identified with the quantum states of the environment. 

In general, the D-wave subspace of reduced density matrix elements (p2)fv involving bilinear 
terms SD*, PD* and DD* is not a decoherence free subspace. Elsewhere we shall show that 
the complete system of SD*, PD* and DD* matrix elements measured on transversely polarized 
target is solvable and yields a unique solution for the D-wave co-evolution amplitudes. However, 
in certain kinematical regions some of these density matrix elements can vanish or be negligible. 



Such is the case of the CERN analysis [42]]. In the kinematic region for dipion mass 1100-1400 
MeV and small momentun transfers t the measurements indicate vanishing of the density matrix 
elements with dipion helicities A = ±2. In this case the solvability of the complete D-wave system 
requires that the co-evolution amplitudes are equal to their average values independent of t. At 
the same time the P-wave amplitudes are small and also nearly equal to their average values. As a 
consequence the same must be true for the S-wave amplitudes, the differences between them being 
three times the diference between P-wave amplitudes (6o| . It is now possible for a trio of helicity 
zero amplitudes S, P° and D° to satisfy approximately the phase and cosine conditions (9.1) and 
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(9.2). Hence we can write 



(ReR°(P T ))H = Re(P?D° T *) ~ |P°||L>?*| cos($(P r °L>°*)) (10.8) 

(ReR°(P T ))<j>l = Re(S T P?*) ~ |£ T ||J**| cos($(S' T P T *)) 

(PeP°(Pr))oo = Re{S T D° T *) ~ |5 r ||£>°*| cos($(,S TJ D°*)) 
The correlation phases were measured in amplitude analysis of data at small momentum transfers 



0.01 < \t\ < 0.20 (GeV/c) and presented at dipion masses above 1180 MeV in Ref. [42|. The phase 
gaps and cosine gaps calculated for these phases are shown in Figure 4. The results are very close 
to and 1, respectively, for both transversities and indicate that the corresponding density matrix 
elements in (10.8) all belong to approximate decoherence free subspace. 

Measurements of the correlation phases determine which density matrix elements belong to 
decoherence or decoherence free subspaces. Both spaces provide information about the nature of 
the interaction of pion creation processes with the environment. 



XI. QUANTUM STATES OF THE ENVIRONMENT. 

The obvious question now is what are the quantum states \e\ > of the environment and how 
do they ensure the diagonal form of co-evolution transversity amplitudes A'jj 1 (£m) = A^ 1 (t)S^ m . 
In a sequel paper [60J] we report on high resolution amplitude analyses of CERN measurements 
of ir~p — ► 7r~7r + n, Tr + n — ► Tr + ir~p and K + n — > K + ir~p on transversely polarized targets below 
dimeson masses ~ 1000 MeV where S- and P-waves dominate. The measured density matrix 
elements do not form a solvable set of equations for the moduli and cosines of correlation phases. 
When the data are supplemented by the assumption that the measured phases satisfy the phase 
relations (9.1), the resulting cosine relations (9.2) provide the missing equations and render the 
system solvable analytically. This amouts to an assumption that the density matrix elements 
describing the S- and P-waves subsystem form a part of the decoherence free subspace. 

In all processes we find two distinct physical solutions for S- and P-wave transversity amplitudes 
A u (i) and A d (j), i,j = 1,2, A = S,P°,P U ,P N leading to 4 distinct final states p f (fl, P, ij). We put 
forward a hypothesis that all four solutions pf(Tt,P,ij) are physical states and that the measured 
final state Pf(£l, P) is a mixed state of these solutions 

2 

p f (n,p)= J2pijPf&,P,ij) (n.i) 

2 

where the probabilities Yl Pij = !■ 

In our next step we associate with the two solutions for transversity amplitudes A u (i) and 
Ad(j), hj = 1,2 two qubit states \i > and \ j >, respectively. Then the hypothesis (11.1) allows 
us to identify the four interacting degrees of freedom of the environment \eg > allowed by the 
condition dimP^(P) < dim Hi (S) dim Hf(S) = (2s p + l)(2s n + 1) = 4 with the four two-qubit 
states \ee >= \i > \j >. Since the transversity amplitudes can possess only one solution at a time, 
we get the diagonal form (8.20) for co-evolution amplitudes A'jj 1 (£m) for any J and A 

A J ^ T {£m) = A J £(££)5 im = A J x ^(ij, ij)8 ij>Vj , = A^(y)%^ (11.2) 
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where 



Kldiv) =< JX V^n\ < ij\U\ij > \0d >= A{ 4 (j) 

Next we identify the probabilities pij in (11.1) with the diagonal terms p mm = Pijij in (8-2). Then 
the equations (8.22) for the moduli have a modified form 

The equations (8.23) for the measured bilinear terms of transversity amplitudes read 

mKXP) - co S (<j>(«t';)) = e *iM*Z!MBt$m (n-s) 



while the equation (8.19) for the measured density matrix elements takes the form 



{ R k)\\> = E 



(11.6) 



In a sequel paper 6l[ we will show that the hypothesis (11.1) is testable experimentally and that 
the probabilities can be measured. 



XII. TIME IRREVERSIBILITY AND VIOLATION OF CPT SYMMETRY 
IN PION CREATION PROCESSES. 

The interacting hadrons are matter particles that carry energy-momentum and spin, and so do 
the two-pion states. The quantum states of the environment 

Pi {E) = Pij,i'j'\i > \i >< i'\ < j'\ (12.1) 

can be thought of as particles carrying quantum entanglement. The unitary co-evolution U of 
the pion creation process with the environment can be interpreted as the scattering of the initial 
hadron state Pi(S) = pi(ir~p, P) with the quantum state Pi{E) of the environment 

Pi(S) ® fH{E) -> u Pf (S,E) (12.2) 

where pf(S,E) is the entangled state of the final hadron system 7r~7r + n with the environment 
E. When observers perform measurements of the final state pf(S,E) they make use of the 
amplitude analysis as an integral part of their measurement and conclude that the observed 
state pf(ir~Tr + n, P) is a mixed state even for pure initial states n~p. This observation informs 
them about the existence of the quantum environment. The observers cannot use their classiclal 
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instruments to observe the environment. It is the pion creation processes that act as non-classical 
instruments which allow observers to access information about the nature of the environment and 
its interactions with hadron processes. 



As we show in a sequel paper [60], the observed state is a mixed state of solutions produced by 
amplitude analysis of data on polarized targets 

2 

p f (Ti~Tr + n,P) = ^ pijPf(7r~Tr + n,P,ij) (12.3) 

The interaction with the environment leads to a "level splitting" of the probability amplitudes 
into solutions A^if) and A^jj), i,j = 1,2 . Such splitting cannot be prepared by the 
observers in the initial state of the time reversed process 7r _ 7r + n — ► ir~p and the pion creation 
process ir~p — > 7r~7r + n is thus time-irreversible. Also, it is not possible for observers to pre- 
pare the entangled state pf(S,E) as an initial state that would undo the entanglement of the 
pion creation process S and the environment E and recreate the initial separable state pi(S)®pi(E). 



According to CPT Theorem, in any local and Lorentz invariant field theory in Minkovski 
spacetime the vacuum expectation values of time-ordered field operators are CPT invariant 



Assuming in addition the unitary S-matrix (1.2), the CPT invariance is extended to the 
observable S'-matrix amplitudes [59']. The CPT invariance then means that the apmlitude A 
which describes a process a + b — > c + d also describes a CPT conjugate process c + d — > a + b. 
This is possible only when masses, widths, spins and isospins of particles and antiparticles are equal. 



Interactions that are invariant under CPT symmetry lead to observables that are invariant 
under the CPT as well. Such is the case for interactions that give rise to hadron masses, widths, 
spins and isospins. Interactions that violate CPT symmetry lead to observables that violate CPT 
symmetry. Interactions (12.2) of pion creation processes with environment violate CPT symmetry 
because of the non-local character of the quantum states Pi(E) and because these states do not have 
well defined antiparticle states. Consequently, the amplitudes A'jj 1 (ij) in (11.3) for the process 

\ir~p > +\i > \j >^ \-K~ir + n > +\i > \j > (12.4) 

do not describe the process 

|7r + 7r~n > +\i > \j >^ \ir + p > +\i > \j > (12.5) 

Since the co-evolution operator U is unitary, the processes (12.2) must be logically time-reversible 
even though they violate CPT symmetry and are time-irreversible from the point of view of 
observers of the reduced states pf(S) who cannot prepare the unknown states pf(S,E). 



It is important to recognize that CPT violating interactions with environment need not 
contradict or affect other CPT invariant interactions and their observables in the same process. 
Thus we may expect the masses and widths of ir~ and ir + to be the same. The CPT violating 
interactions will manifest themselves in other observable aspects of the pion creation processes. In 
a series of sequel papers on high resolution amplitude analyses we will show that the interaction 
with environment manifests itself in mixing of scalar and vector resonances below ~ 1000 MeV in 
both ttN — > tc~tt + N [60] and K + n — >• K + n~p processes, and in violation of the Generalized Bose- 
Einstein symmetry in ir~p — > -ir~ir + n [6o| . In previous low resolution analyses of irN — > it~it + N 
the presence of p°(770) resonance in the S*-wave amplitudes was misinterpreted as evidence for a 
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scalar resonance cr(770) 



U EH, 0, 0,(54]. 



There is no energy-momentum of the quantum state Pi(E) of the environment involved in the 
interaction with pion creation processes. The interacting hadrons conserve their energy-momentum 
and there is no exchange of energy-momentum with the environment, in agreement with the original 
proposal by Hawking for particle processes interacting with quantum fluctuations of space-time 
metric [8j, |9j. Instead, the interaction with environment is a non-dissipative (dephasing) process 
and involves change of quantum entanglement of quantum states of the environment and the 
produced hadrons. The entanglement changing interaction modifies the entanglement content of 
the quantum state of the environment to Pf(E) = Trs(pf(S, E)) and stores quantum information in 
dimeson states observable as mixing of scalar- vector resonances [60] and as violation of Generalized 
Bose-Einstein symmetry 60] . This last observation seems to confirm the change of entanglement of 
particle-antiparticle pairs called "w-effect" and recently predicted by Bernabeu, Mavromatos and 
Sarkar to arise in CPT violating interactions of maximally entangled particle-antiparticle pairs 
with space-time foam 2ll. |22|. 



XIII. CONCLUSIONS. 



It is the task of spin physics to use known initial spin states pi to measure and to investigate 
the final spin states pf to learn about new aspects of the dynamics of hadron interactions. It is 
therefore not surprising that the CERN measurements of pion creation processes on polarized 
targets provide information about the validity of the unitarity assumptions (1.1) and (1.2). The 
observed violation of the unitarity conditions (4.15) and (4.16) means that pure initial states can 
evolve into mixed final states in these processes. This leads to a new view of pion creation processes 
as open quantum systems interacting with a quantum environment in a unitary co-evolution. 
The non-unitary evolution from the initial state to the observed final state is described by 
Kraus representation that leaves invariant the form of measured angular distributions in terms of 
environment averaged density matrix elements. 

The measured process is time- irreversible and violates CPT symmetry. Due to the non-local 
character of the interaction with the environment and due to the fact that the quantum states 
Pi{E) of the environment do not have well defined antiparticle states, the co-evolution process itself 
violates CPT symmetry. The interaction with the environment manifests itself in scalar- vector 
mixing and in violation of Bose-Einstein symmetry. 

The CERN measurements of pion creation processes on polarized targets opened a window on a 
entirely new class of phenomena which result from their interactions with a quantum environment. 
Following Hawking's proposals, these phenomena can be understood as low energy manifestations 
of quantum gravity. In this view, the observed quantum environment originates in quantum gravity 
and the pion creation processes act as non-classical instruments that make quantum gravity an 
observable reality. 
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APPENDIX A: LUTZ-RYBICKI TABLES OF DENSITY MATRIX ELEMENTS. 



In their 1978 paper [32], Lutz and Rybicki tabulated expressions for all 16 kinds of angular 
density matrix elements {R 3 k )\{i , 3 = 0, 1, 2, 3, k = u, x, y, z in terms of helicity amplitudes, helic- 
ity amplitudes with definite naturality, nucleon transversity amplitudes and nucleon transversity 
amplitudes with definite naturality. Their work has been used in all amplitude analyses of CERN 
data on ttN — ► tt~tt + N and K + n — > K + ir~p on polarized targets 41, 43, 44, 48, 4^, 51, 5^] but was 



not published. Since their paper is no longer available, we reproduce their tables for helicity am- 
plitudes and nucleom transversity amplitudes with definite naturality, which are used in this work, 
in Figures 5 and 6, respectively. In our notation their nucleon helicity nonflip and flip amplitudes 
read 

NL = H( +fi+ , F m = H J x+fi _ (Al) 
The relations for nucleon transversity amplitudes with definite naturality read 

^gi n = U( = Ui iUl ^hi n = U J x = Ui d (A2) 

N gl = N{ = N( 4 , N hl n = N( = N{ u 

In their notation, the parity relations (2.52) for transversity amplitudes with definite naturality 
are 

u gl m = +(-l) mC/ ^, U hl m = +(-l) ml V m (A3) 

9-m = ~ (~ -0 9mi h—m = ~ (~ -Q 

They use lower case jf and mm' while we use J J' and AA' with the same meaning. 

Nucleon helicity amplitudes with definite t-channel naturality defined in (2.44) and (2.45) are 
better suited for certain theoretical studies. Their characteristic feature is pion exchange dominance 
of unnatural helicity flip amplitudes at small momentum transfers t. In the companion paper [§3] 
we show how conversion of transversity amplitudes into helicity amplitudes with definite naturality 
determines the relative phase of transversity amplitudes with opposite nucleon transversity and thus 
recoil nucleon polarization. In Figure 7 we reproduce Lutz-Rybicki Table lb for nucleon helicity 
amplitudes with definite t-channel naturality. In our notation their amplitudes read 

U < = t^ +) o+i N < = N X+,o+ (A4) 

U fj _ TjJ N fj _ AT J 

Their transversity amplitudes (A2) are related to their helicity amplitudes (A4) by relations 

9m = "^2^ n ^ m * frn)y ^m = ~~^J^ n ™ ~ * f™>) (A5) 

N 3 = J_(N j _ -iV f j x N h j = J_(N j , aN f j x 
\/2 JmJi n m 

V = N K = o 

These are the same relations as (2.49) and (2.50). To achieve this conformity was the reason why 
we omitted the factor i in (2.47). 
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FIG. 5: Angular density matrix elements expressed in terms of helicity amplitudes. The spin indices jf 
which always go with mm' have been omitted in the amplitudes. Reproduced Table la of Ref. [32 1. 
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FIG. 6: Angular density matrix elements expressed in terms of nucleon transversity amplitudes with definite 
i-channel naturality. The spin indices jj' which always go with mm' have been omitted in the amplitudes. 
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FIG. 7: Angular density matrix elements expressed in terms of nucleon helicity amplitudes with definite 
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